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Abstract 


We demonstrate a new computational illumination technique that achieves large 
space-bandwidth-time product, for quantitative phase imaging of unstained live sam¬ 
ples in vitro. Microscope lenses can have either large field of view (FOV) or high 
resolution, not both. Fourier ptychographic microscopy (FPM) is a new compu¬ 
tational imaging technique that circumvents this limit by fusing information from 
multiple images taken with different illumination angles. The result is a gigapixcl- 
scale image having both wide FOV and high resolution, i.e. large space-bandwidth 
product (SBP). FPM has enormous potential for revolutionizing microscopy and has 
already found application in digital pathology. However, it suffers from long acquisi¬ 
tion times (on the order of minutes), limiting throughput. Faster capture times would 
not only improve imaging speed, but also allow studies of live samples, where motion 
artifacts degrade results. In contrast to fixed (e.g. pathology) slides, live samples 
are continuously evolving at various spatial and temporal scales. Here, we present a 
new source coding scheme, along with real-time hardware control, to achieve 0.8 NA 
resolution across a 4x FOV with sub-second capture times. We propose an improved 
algorithm and new initialization scheme, which allow robust phase reconstruction over 
long time-lapse experiments. We present the first FPM results for both growing and 
confluent in vitro cell cultures, capturing videos of subcellular dynamical phenomena 
in popular cell lines undergoing division and migration. Our method opens up FPM 
to applications with live samples, for observing rare events in both space and time. 



1 Introduction 


In vitro microscopy is crucial for studying physiological phenomena in cells. For 
many applications, such as drug discovery [l], cancer cell biology |2| and stem cell 
research |3|, the goal is to identify and isolate events of interest. Often, these events 
are rare, so automated high throughput imaging is needed in order to provide sta¬ 
tistically and biologically meaningful analysis |4|. Thus, an ideal technique should 
be able to image and analyze thousands of cells simultaneously across a wide field 
of view (FOV). To observe dynamical processes across various spatial and temporal 
scales, both high spatial and high temporal resolution are required. Existing high 
throughput imaging techniques [5]-[l0 cannot meet the space-bandwidth-time prod¬ 
uct (SBP-T) required for widefield in vitro applications. Here, we introduce a new 
computational microscopy technique with both high spatial and temporal resolution 
over a wide FOV. 

Our method is an extension of Fourier ptychographic microscopy (FPM) |7|, which 
overcomes the physical space-bandwidth product (SBP) limit. Instead of choosing be¬ 
tween large FOV or high resolution, FPM achieves both by trading acquisition speed. 
Illumination angles are scanned sequentially with a programmable LED array source 
(Fig-0. while taking an image at each angle. Tilted illumination samples different 
regions of Fourier space, as in synthetic aperture 11,12 and structured iliumina¬ 
imaging. Though the spatial resolution of each measurement is low, the 
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tion 

images collected with high illumination angles (dark field) contain sub-resolution in¬ 
formation. The reconstructed image achieves resolution beyond the diffraction limit 
of the objective - the sum of the objective and illumination NAs. Distinct from syn¬ 
thetic aperture, FPM does not measure phase directly at each angle, but instead uses 
nonlinear optimization algorithms (7,15 17 similar to translational diversity 18 19 


and ptychography 20,21 . Conveniently, the LED array coded source used here is 


implemented as a simple and inexpensive hardware modification to an existing mi¬ 
croscope. Our proposed method uses efficient source coding schemes to reduce the 
capture time by several orders of magnitude. 

The standard approach for large SBP microscopy is slide scanning, in which the 
sample is mechanically moved around while imaging with a high resolution objective 
to build up a large FOV. Scanning limits throughput 1221 and is unsuitable for in vitro 
imaging of dynamic events. Instead of starting with high resolution and stitching 
together a large FOV, FPM starts with a large FOV and stitches together images 
to recover high resolution. A major advantage is its ability to capture the required 
set of images with no moving parts, by varying illumination angles. In addition, 
scanners use high magnification objectives with short depth of field (DOF), thus 
requiring extensive auto-focusing mechanisms. In contrast, FPM provides robustness 
to focus errors because its DOF is longer than that which would be provided by a 
high-magnification objective of equivalent NA 17 . Further, out of focus images can 
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Figure 1: Source-coded Fourier ptychographic microscopy (FPM) captures large 
space-bandwidth product (SBP) images in under 1 second. (A) Experimental setup 
is a microscope with an LED array source and a wide FOV 4x (0.2 NA) objective. 
Multiple images are captured with coded illumination in order to reconstruct higher 
resolution (up to 0.8 NA). (B) Comparison of illumination schemes in terms of space, 
bandwidth and acquisition time. Sequential FPM scans through each LED, achieving 
large SBP at a cost of speed. Our source-coded FPM implements hybrid patterning 
to achieve the same SBP with sub-second acquisition time. (C) The number of images 
required for source-coded FPM (blue) grows more than 8 x slower than for sequential 
FPM (red) as the final resolution increases. (Solid lines: theoretical, points: our LED 
array). 


be digitally refocused |7|. 

The main limitation of FPM for real-time in vitro applications is long acquisition 
times and large datasets. Not only are a large number of images (~200) captured 
for each reconstruction, but also long exposure times are needed for the dark field 
images. To shorten exposure times, we first built a custom hardware setup consisting 
of high-brightness LEDs and fast control circuitry (see Methods 4.1). This allows us to 
reduce total acquisition time by ~ 50 x, enabling a full 173 LED scan to capture 0.96 
gigapixels of data in under 7 seconds. While this speed is already suitable for many in 
vitro applications (e.g. cell division processes), fast subcellular dynamics (e.g. vesicle 
tracking) require sub-second capture to avoid severe motion blur. Since the camera’s 
data transfer rate is the limiting factor, large SBP images cannot be captured in 
less than 1 second unless we can reduce the data requirements. Furthermore, for 
studies requiring time-series measurements, the large data requirement (~1 gigapixel 
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per dataset) poses severe burdens on both storage and processing. 

In order to reconstruct large SBP from fewer images, we need to fundamentally 
change our capture strategy. To do this, we eliminate redundancy in the data by 
designing new source coding schemes. The redundancy arises because of a 60% over¬ 
lap requirement in Fourier space for neighboring LEDs |T, 23]. This means that we 
must capture ~10x more data than we reconstruct if using ‘sequential’ FPM. Our 
angle-multiplexing scheme instead turns on multiple LEDs simultaneously for each 
measurement, allowing better coverage of Fourier space with each image. Without 
eliminating the overlap, we therefore £11 in Fourier space faster. Previous work em¬ 


ployed a random coding strategy across both brightfield and dark field regions 16 


One problem with this scheme is that the brightfield images and dark field images 
have large differences in intensity (~10-100x). Considering Poisson noise statistics, 
this means that images with mixed brightfield and dark field components may suffer 
from the dark field signal being overwhelmed by the brightfield noise. As a result, 
we significantly improve our multiplexing results here by separating brightfield from 
dark field LEDs. Furthermore, we note that asymmetric illumination-based differ¬ 
ential phase contrast (DPC) 124, 25 provides a means for recovering quantitative 


phase and intensity images out to 2x the objective NA with only 4 images. Hence, 
there is no need to individually scan the brightfield LEDs. Our new method, termed 
source-coded FPM, uses a hybrid illumination scheme: it first captures 4 DPC images 
(top, bottom, left, right half-circles) to cover the brightfield LEDs, then uses random 
multiplexing with 8 LEDs to fill in the dark field Fourier space region (Fig. [jj3). 

Source-coded FPM approaches the theoretical limit for SBP-T [26,2T|, set by the 
camera’s data transfer rate. We achieve SBP-T of 46 megapixels/second, calculated 
according to [28j (4x FOV and 0.8 NA captured in 0.8 seconds). Though the data rate 
of our camera is larger (138 megapixels/second), we leave some redundancy in order to 
ensure robust algorithm convergence. Reconstruction algorithms that explicitly take a 
priori information into account, such as sparsity based methods 29 30 , could further 
improve capture speed; here, for generality, we choose not to make any assumptions 
on the sample. Both sequential FPM and source-coded FPM require the number of 
images in the dataset to grow quadratically with improved final resolution. This is 
due to the coverage area in Fourier space increasing in proportion to the square of the 
final NA. However, our source-coded method has a slower growth rate (Fig. R) and 
stays flat for less than 2x resolution improvement. Conveniently, these techniques can 
flexibly trade off FOV, resolution and acquisition time by choosing the illumination 
angle range. 

For in vitro applications, stain-free and label-free is particularly attractive be¬ 
cause it is non-invasive and non-toxic. FPM provides both intensity and phase 
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which contain morphological and cell mass 1321 information that can be used for 


quantitative study 33,34 . In order to achieve accurate, high-quality results for un¬ 


stained samples, we needed to make several modifications to the FPM algorithm. 
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It is well known that non-convex problems such as phase retrieval will often get 
stuck in local minima 35-37 ; the best way to avoid this is to provide a good initial 


guess (8,36,38 . Previous work in FPM uses a low-resolution intensity image as the 
initial guess (7 15 16,31 . Stained samples with strong intensity variations thus recon¬ 
struct successfully, since the intensity-only initialization is close to the actual solution. 
However, unstained samples are phase objects and so the intensity-only initial guess 
does not provide a good starting point. Furthermore, FPM does not measure phase 
directly at each angle, and the phase contrast provided by the asymmetric illumina¬ 
tion results in uneven sensitivity to phase at different spatial frequencies. A detailed 
phase transfer function analysis 25 based on the weak-object (Born) approxima¬ 


tion 39,40 shows that low-frequency phase information is poorly captured, since it 
results only from illumination angles that are close to the objective NA. Thus, low 
spatial frequency phase information is more difficult to reconstruct than high spatial 
frequency phase information, contrary to the situation for intensity reconstructions. 
To improve our reconstruction, we use a linearly approximated phase solution based 
on DPC deconvolution 25 as a close initial guess for spatial frequencies within the 


2 NA bandwidth. We then run a nonlinear optimization algorithm to solve the full 
phase problem (see Methods 4.2), resulting in high-quality phase reconstructions with 
high resolution (Fig. (2}A) and good low spatial frequency phase recovery (Fig. (2^3) . 

We demonstrate our new source-coded FPM by reconstructing large SBP videos 
of popular cell types in vitro on a petri dish, for both growing and confluent samples. 
We observe sub-cellular and collective cell dynamics happening on different spatial 
and temporal scales, allowing us to observe rare events in both space and time, due 
to the large space-bandwidth-fzme product and flexible tradeoff of time and SBP. 


2 Results 

2.1 Validation with stained and unstained samples 

A major advantage of quantitative phase imaging is that it can visualize transparent 
samples in a label-free, non-invasive way. Figure U4 compares our reconstructions 
before and after staining, for the same fixed human osteosarcoma epithelial (U20S) 
sample. With staining, the intensity image clearly displays detailed subcellular fea¬ 
tures. Stained samples also display strong phase effects, proportional to the local 
shape and density of the sample. Without staining, the intensity image contains very 
little contrast and is nearly invisible; however, the phase result clearly captures the 
sub-cellular features. Due to the strong similarity between stained intensity and un¬ 
stained phase, it follows that quantitative phase may provide a valid alternative to 
staining. 

To demonstrate the importance of using a good initial guess to initialize the phase 
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Figure 2: Large SBP reconstructions of quantitative phase and intensity. (A) Phase 
reconstruction across the full FOV of a 4x objective with 0.7 NA resolution (sam¬ 
ple: U20S). A zoom-in is shown at right, with comparison to reconstructions of the 
same sample before and after staining. (B) Our improved FPM algorithm provides 
better reconstruction of low-frequency phase information. A zoom-in region shows 
comparisons between phase reconstructions with and without our differential phase 
contrast (DPC) initialization scheme. (C) To validate our source-coded FPM results, 
we compare to images captured with a 40 x objective having high resolution (0.65 
NA) but small FOV (sample: MCF10A), as well as sequential FPM. (D) We simulate 
a phase contrast image and compare to one captured by a high-resolution objective 
(0.65 NA, 40x). 

recovery for unstained samples, we compare the FPM results both with and without 
our DPC initialization scheme (Fig. §3). Both achieve the same 0.7 NA resolution, 
with high spatial frequency features (e.g. nucleus and hlopodia) being clearly recon¬ 
structed, as expected. However, without DPC initialization, low frequency compo¬ 
nents of the phase are not well recovered, resulting in a high-pass filtering effect on 
the reconstructed phase, much like Zernike phase contrast. With DPC initialization, 
low frequencies, which describe the overall height and shape of the cells, are correctly 
recovered. 

Next, we verify the accuracy of the recovered phase values by comparing to images 
captured directly with a higher resolution objective (40x 0.65 NA). The reconstruc¬ 
tion resolution for our method matches its expected value (0.7 NA), as seen in Fig. [2]C, 
which shows images of fixed human mammary epithelial (MCF10A) cells. To validate 
our quantitative phase result, we compare to that recovered from a through-focus in¬ 
tensity stack captured with the high-resolution objective, then input into a phase 
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Figure 3: Time-lapse large SBP phase reconstruction of unstained HeLa cells un¬ 
dergoing division. (A) Sample raw data and Fourier coverage using sequential FPM 
(173 images), with acquisition time of 7 seconds per frame. (B) One frame of the 
full FOV phase reconstruction using a 4x objective and achieving 0.8 NA resolution. 
(C) Several frames of reconstructed video (see Video 1) from a zoom-in of one small 
area of confluent cells in which one cell is dividing into multiple cells. (D) Automated 
cell segmentation result for the full FOV phase image, with ~3,400 cells successfully 
identified. (E) Calculated dry mass for each of the labeled cells in the zoom-in re¬ 
gion over 4 hours at 2 minutes intervals. At right is a histogram of the background 
fluctuations in an area with no cells. 


retrieval algorithm |4l] (Fig. [2]C and Methods 4.3). Quantitative phase can further 
be used to simulate popular phase contrast modes such as DIC and Zernike phase 
contrast (PhC) |42|. In Fig. 03. we compare an actual PhC image to that simulated 
from our method’s reconstructed result (see Methods 4.4). Since the PhC objec¬ 
tive uses annular illumination (Ph2, 0.25 NA) to achieve resolution corresponding to 
0.9 NA, it should have slightly better resolution than our reconstruction. Note that 
the simulated PhC image is effectively high-pass filtered phase information, so small 
details appear with better contrast. 


2.2 Fast sequential FPM video of HeLa cells dividing in vitro 

Figure [3] shows a few frames from a time-lapse video of human cervical adenocar¬ 
cinoma epithelial (HeLa) cell division process over the course of 4 hours, as well as 
an automated quantitative analysis of cell dry mass evolution. This data was cap¬ 
tured using our improved sequential FPM (173 images) - sample raw images and a 
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schematic of the Fourier space coverage are shown in Fig. [3j\ . One time frame of the 
full FOV phase reconstruction is shown in Fig. [3j3. and a few frames of the video 
for a single zoom-in are shown in Fig. [3p. In this region, one cell is undergoing mi¬ 
tosis and dividing into 4 cells, during which the cells detach from the surface and 
become more globular (see Video 1). This situation, if imaged with a high magnifi¬ 
cation objective, would result in the detached cells moving out of the focal plane and 
blurring. However, FPM provides a longer DOF than a high magnification objective 
with equivalent NA, so the entire sample stays in focus across the FOV. Subcellular 
features are visible and the dynamics of actin fila m ent formation can be tracked over 
time. We show only phase results (omitting intensity), since samples are unstained 
and have little intensity contrast. 

Since our reconstructed video contains ~20 gigapixels of quantitative phase data 
with ~3,000 cells in each frame, the practical extraction of relevant information re¬ 
quires automated analysis. It is well known that PhC/DIC images cannot be auto¬ 
matically segmented (due to the lack of low spatial frequency information); however, 
our quantitative phase results do not have this problem. Using an automated cell seg¬ 
mentation software (Cellprohler |43|) applied directly to the full FOV phase result, we 
first segment each frame of the video to find each of the ~3,400 cells (Fig. |3]D). Next, 
we compute each cell’s dry mass 32,34 over time through the division process, with 
a few sample cells plotted in Fig. (3F. Dry mass is calculated by integrating phase over 


each segmented cell region for each time frame (see Methods 4.5). Note that the au¬ 
tomated cell segmentation and dry mass calculation are sensitive to the quality of the 
phase result, and often fail when low spatial frequencies are not well reconstructed; 
hence, the DPC initialization and other algorithm improvements implemented in this 
work are crucial for automated quantitative studies. 


2.3 Source-coded FPM video of neural stem cells in vitro 

Source-coded FPM can observe samples across long time scales (up to 4.5 hours) with 
sub-second acquisition speed (1.25 Hz). An example frame from a reconstructed large 
SBP phase video of adult rat neural stem cells (NSCs) is shown in Fig. [4jB, in a petri 
dish. We use source-coded FPM to achieve the same SBP as in sequential FPM (0.8 
NA resolution across a 4x FOV), but with only 21 images (sample raw images shown 
in Fig. [4]W) , as opposed to 173. As a result, we significantly decrease the capture time, 
from 7 seconds to 0.8 seconds. This alleviates motion artifacts that would otherwise 
blur the result due to sub-cellular dynamics that happen on timescales shorter than 
7 seconds. Since our source-coded FPM reduces the number of images needed, we 
can also capture longer video sequences without incurring data management issues. 
Hence, we create videos of both fast-scale dynamics and slow-scale evolution of NSCs 
with details both on the sub-cellular level and across the entire cell population. For 
example, while vesicle transport and other organelle motions tend to occur on a short 


7 







/l. c*.• *• * 200|lim 


' ; t *. ^ . „r. 3 &. v 

*■ . ^ *, . 4 * a . 

• • , ■*. . . \ . 
v ’ * 1 *s~ **. •’ % ‘ t , . « . 

7 •• > •'.♦.* • .. .7. / 

V.'- ’ vv.; v''■'■'( '■ •;• > 


4 DPC 17 multiplexed 
images dark field images) I 

Fourier space + K 


& 


objective’s 

bandwidth 


•• «> . f v 

. \ - ' > ; - ~ ' • A 
- * 

v/ . i* -v *; * 


~f *'*'*•'* ‘ * : ! * 




• ■ »/, 



I c 



Figure 4: Large SBP phase video reconstructions for observing multi-scale temporal 
dynamics of in vitro NSCs with high SBP and acquisition time of 0.8 seconds per 
frame. (A) Our source-coded FPM captures 4 brightfield images and 17 multiplexed 
dark held images. (B) Full FOV phase reconstruction using a 4x objective and 
achieving 0.8 NA resolution. (C) Sample frames of reconstructed video (see Video 2) 
for a zoom-in of one small area. Top: successive frames at the maximum frame rate 
(1.25 Hz), Bottom: sample frames across the longer time-lapse (4.5 hours at 1 minute 
intervals). 

time scale and at a small length scale (see Video 2), NSC differentiation into neuronal 
and glial lineages occurs on a longer time scale (he. days) and at a larger length 
scale. In the middle time scale, we observe NSCs sensing their micro-environment by 
retracting and extending processes, reorganizing their cytoskeletons, migrating, and 
maturing their axonal and dendritic processes. 
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Figure 5: Motion blur degrades effective resolution in live dynamic samples. Recon¬ 
structed phase of live samples using different capture schemes with the same nom¬ 
inal spatial resolution (0.8 NA), but different acquisition times. As capture speed 
increases, more details about subccllular dynamics become visible due to reduced 
motion blur. Two fast dynamical processes in MCF10A cells: (A) subcellular fiber 
motion and (B) vesicle transport (Video 3) are blurred out when acquisition times are 
longer than 1 second. Our source-coded FPM achieves sub-second capture, revealing 
more details, yet not as clearly as DPC, which has the fastest capture time. (C) 
and (D) results for NSCs, which exhibit slower dynamics than the MCF10A cells. 
Sequential FPM blurs out most subcellular features; however, our source-coded FPM 
is able to capture details without motion artifacts (Video 4). 

3 Discussion 

With live samples imaged in vitro, dynamics create motion blur artifacts that can 
destroy the resolution improvements gained by large SBP methods. Thus, the final 
effective resolution is always coupled with acquisition speed and sample-dependent 
motion. In general, smaller subccllular features tend to move at faster speed; hence, 
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we find that capture times on the order of minutes always incur motion blur. For 
this reason, faster acquisition is essential for in vitro applications with typical cell 
types. To demonstrate this point, we compare the results of our source-coded FPM, 
sequential FPM with and without real-time hardware controls, and high magnifica¬ 
tion DPC. The DPC method is considered our benchmark, since it achieves faster 
acquisition speeds, avoiding motion blur, albeit with a small FOV. In Fig. [5j we show 
zoomed-in phase reconstructions for two cell types with each capture scheme, from 
slowest to fastest. In each case, the final result has a nominal NA of 0.8 and so each 
of these results should have the same resolution. However, it is clear that the slower 
capture schemes blur out features, rendering much of the small-scale structure and 
dynamics invisible. 

The MCF10A cells in Fig. [5)4, B exhibit the fastest dynamics observed, due to 
rapid shuttling of small vesicles and fluctuations of cytoskeletal fibers (e.g. actin 
filaments, microtubules). Using sequential FPM with 60 second capture time, we 
retain almost no details about the structure of the fibers (Fig. (5)4). Even with our 
faster hardware setup, requiring only 7 seconds acquisition time, the fibers are still 
completely blurred out. By switching to our source-coded FPM, we obtain the same 
SBP as the previous two schemes, but with only 0.8 seconds acquisition time. Now, 
fiber cytoskeletal details become more discernible, though there is still some motion 
blur, as compared to the DPC result (Video 3). Hence, even our source-coded FPM is 
missing some information in this case. In Fig. [5j3 we zoom in on another fast process 
of vesicle transport, where we find that our method captures more of the dynamics 
than in the fiber fluctuation case, still with some blurring. In Video 3, we show time- 
lapse videos over 8 hours to observe cell migration and proliferation with subcellular 
details, captured using our source-coded FPM. Thus, the trade-off in FOV and time 
may be more appropriate, depending on what one aims to observe. 

In Fig. [5p,D, example results for two different NSCs are shown. In both cases, 
sequential FPM results in significant motion blur, particularly along thin, extended 
processes and within intracellular vesicles and organelles. In contrast, source-coded 
FPM is able to accurately capture the full details of the sample without motion blur 
artifacts. While it is difficult to compare directly because the live cells were moving in 
between capture schemes, in general we can say that most vesicle transport, retraction 
and extension of processes and other organelle motion can be clearly visualized using 
source-coded FPM (see Video 4). 

The flexibility of our system in trading off FOV, resolution and time means that 
experiments can be tailored to the sample. For example, the HeLa and NSC samples 
shown here display slower sub-cellular dynamics, so are more suited to high SBP 
imaging with our current scheme, whereas the MCF10A sample necessitates a trade¬ 
off of SBP in order to capture data with sufficient speed. When dynamics are on 
timescales faster than 0.8 seconds, one should reduce the number of images captured, 
either by sacrificing FOV (using a larger NA objective) or by reducing the resolution 
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improvement factor (using a smaller range of LEDs). In the limit, one can eliminate all 
the dark held LED images and simply implement DPC for maximum speed of capture. 
However, for most biological dynamics that are studied in vitro (e.g. differentiation, 
division, apoptosis), we End that sub-second acquisition is sufficient, and so the SBP 
should be maximized within this constraint. 

In summary, we have demonstrated a high-speed, large SBP microscopy technique, 
providing label-free quantitative phase and intensity information. Our source-coded 
FPM method overcomes the limitations of existing large SBP methods, permitting 
fast, motion-free imaging of unstained live samples. This work opens up large SBP 
imaging to high throughput in vitro applications across a large range of both spa¬ 
tial and temporal scales. A gallery of interactive full FoV high resolution images 
from our experimental system can be found at http: //www. gigapan. com/profiles/ 
WallerLab_Berkeley. 


4 Methods 

4.1 Experimental setup 

We place a custom-built 32x32 surface-mounted LED array (4mm spacing, central 
wavelength 513nm with 20nm bandwidth) placed ~70mm above the sample (Fig.[lj4), 
replacing the microscope’s standard illumination unit (Nikon TE300). All LEDs 
are driven statically using 64 LED controller chips (MBI5041) to provide indepen¬ 
dent drive channels. A controller unit based on ARM 32-bit Cortex™ M3 CPU 
(STM32F103C8T6) provides the logical control for the LEDs by the I 2 C interface at 
5MHz, with LED pattern transfer time of ~320/is. The camera (PCO.edge 5.5, 6.5/im 
pixel pitch) is synchronized with the LED array by the same controller via two coaxial 
cables which provide the trigger and monitor the exposure status. All raw images are 
captured with 14ms exposure times. We experimentally measure the system frame 
rate to be ~25Hz for capturing full-frame (2560 x 2160) 16 bit images. The data 
are transferred to the computer via CameraLink interface. All in vitro experiments 
are performed in petri dishes placed inside a temperature and CO 2 controlled stage 
mounted incubator (I 11 Vivo Scientific). 


4.2 Large SBP quantitative phase and intensity reconstruc¬ 
tion 


Our new FPM reconstruction algorithm can be described in two steps. First, we cal¬ 
culate a low-resolution initialization based on DPC. Next, we implement our quasi- 
Newton’s method iterative reconstruction procedure 16 to include the higher-order 
scattering and dark held contributions, for 3-5 iterations. In our source-coded FPM, 
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the 4 brightfield images are directly used in the deconvolution-based DPC recon¬ 


struction algorithm 25 to calculate phase within 2x the objective’s NA. The initial 


low-resolution intensity image is calculated by the average of all brightfield images 
corrected by the intensity falloffs |44J and then deconvolved by the absorption trans¬ 
fer function 25 . In onr improved sequential FPM algorithm, the 4 DPC images are 


numerically constructed by taking the sum of single-LED images corresponding to 
the left, right, top, bottom half-circles on the LED array, respectively. 

In the reconstruction, we divided each full FOV raw image (2560x2160 pixels) 
into 6x5 sub-regions (560x560 pixels each), with 160-pixel overlap on each side of 
neighboring sub-regions. Each set of images were then processed by our algorithm 
above to create a high resolution complex-valued reconstruction having both inten¬ 
sity and phase (2800x2800 pixels). Finally, all high resolution reconstructions were 
combined using the alpha-blending stitching method [7j to create the full FOV high 
resolution reconstruction. Using a desktop computer (Intel i7 CPU), the processing 
time for each sub-region was ~30s in Matlab. The total processing time for each full 
FOV was rs_/ 20min. 


4.3 Phase reconstruction from through-focus intensity stack 

We capture intensity images with a high-resolution objective (40x 0.65 NA) while 
moving the sample axially to 17 exponentially spaced positions from -64/im to 64yum 
using a piezostage (MZS500-E-Z, Thorlabs). The images are then used to recon¬ 
struct phase (Fig. [2p) using a transport of intensity type algorithm based on spatial 


frequency domain fitting 41 


4.4 Simulation of phase contrast images 

Our phase contrast simulation (Fig. [2jD) fully accounts for the partially coherent 
annular illumination (inner NA = 0.25, outer NA = 0.28, measured experimentally) 
and apodized phase contrast pupil 45 . The pupil consists of a 7r/2-phase shifting 


ring with 75% attenuation (size matching the source), and two apodization rings 
with 50% attenuation (widths are calculated based on 45|). In our simulation, a 
tilted plane wave from each source point shifts the sample’s spectrum in the Fourier 
space, which is then filtered by the pupil function before calculating the intensity 
image in the real space. The phase contrast image is the incoherent sum of all the 


intensity contributions from all the points on the annular source 46 


4.5 Image segmentation and cell dry mass calculation 


Image segmentation for each frame is performed by Cellprobler 43 open-source soft¬ 


ware, which implements a series of automated operations including thresholding, wa- 
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tershedding and labeling, to return a 2D map containing segmented regions repre¬ 
senting different cells. The 2D maps are then loaded into Matlab to extract the phase 
within each individual cell. The total dry mass for each cell is calculated as the sum 
of the dry mass density 32 . The dry mass density p is directly related to phase 


(j) by p = gwhere A is the center wavelength, 7 = 0.2mL/g is the average of 


reported values for refractive increment of protein |47|. The background fluctuations 
are characterized from a region without any cells (white square in Figure |3]D) and 
having an area similar to the average cell size. Background fluctuations are shown 
in the black curve and histogram in Figure [3]E. We achieve a standard deviation of 
1.5pg, in units of dry mass, indicating good stability of the phase measurement. 


4.6 Sample preparation 

HeLa cells were cultured with DMEM (Dulbeccos Modified Eagles Medium) supple¬ 
mented with 10% FBS (fetal bovine serum), Glutamine and penicillin/streptomycin. 
The cells were plated on a p75 flask and cultured in a 37°C incubator with 5% 
C0 2 • The confluent cells were treated with 0 . 2 % trypsin and passed as 1:8 ratio. 
Two drops of trypisized cells were then added into a poly d-d-lysine coated 35mm 
MatTek glass-bottom plate with 2mL medium. After 6 hours incubation, the cells 
became fully attached to the plate. U20S cells were prepared using the same pro¬ 
cedure. They were fixed in 4% formaldehyde at room temperature for 10 min and 
later stained with 1% Toluidine Blue O (TBO) at room temperature for 5 min and 
rinsed in three changes of double-distilled water. Adult rat NSCs were isolated from 
the hippocampi of 6 -week-old female Fischer 344 rats. To promote adhesion, tissue 
culture polystyrene plates were first coated with 10/rg/mL poly-L-ornithine (Sigma) 
overnight at room temperature, followed by 5/rg/mL of laminin (Invitrogen) overnight 
at 37°C. NSCs were cultured in monolayers in 1:1 DMEM/F12 high-glucose medium 
(Life Technologies), supplemented with N-2 (Life Technologies) and 20ng/mL re¬ 
combinant human FGF-2 (Peprotech). Media was changed every other day, and cells 
were subcultured with Accutase upon reaching ~80% confluency. MCF10A cells were 
cultured in DMEM/F12 (Invitrogen), supplemented with 5% horse serum (Invitro¬ 
gen), 1% penicillin/streptomycin (Invitrogen), 0.5/ig/mL hydrocortisone (Sigma), 100 
ng/mL cholera toxin (Sigma), 10/rg/mL insulin (Sigma) and 20ng/mL recombinant 
human EGF (Peprotech). Media was changed every other day, and cells were passed 
with Trypsin upon reaching ~80% confluency. I 11 preparation for imaging, cells were 
washed once with PBS (phosphate buffer saline) and detached with either Accutase 
or Trypsin. 300,000 cells were seeded onto 35 mm glass-bottom microwcll dishes 
(MatTek) and allowed to attach. 
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